LAEFF-95/004 
LA-UR-95-1055 

Ising Model formulation of Large Scale Dynamics: 
Universality in the Universe. 

J. Perez-Mercader* , T. Goldman''', D. Hochberg*^ and R. Laflamme^ 

*Laboratorio de Astrofisica Espacial y Fisica Fundamental 
Apartado 50727 
28080 Madrid, Spain 
' Theoretical Division 
Los Alamos National Laboratory 
Los Alamos, New Mexico 87545, USA 
ft Departamento de Fisica Teorica 
Universidad Autonoma de Madrid 
Cantoblanco, 28049 Madrid, Spain 



Abstract 



The partition function of a system of galaxies in gravitational interaction 
can be cast in an Ising Model form, and this reformulated via a Hubbard- 
Stratonovich transformation into a three dimensional stochastic and classical 
scalar field theory, whose critical exponents are calculable and known. This 
allows one to compute the galaxy to galaxy correlation function, whose non- 
integer exponent is predicted to be between 1.530 and 1.862, to be compared 
with the phenomenological value of 1.6 to 1.8. 

PACS: Self-Gravitating Systems: 04.40; Cosmology: 98.80; Renormalization, - phase 
transitions: 64. 60 A. 
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An intergalactic tourist admiring the Universe at the largest scales would perceive it 
as something akin to a three dimensional 'salt-and-pepper' pattern which, when projected 
onto a two dimensional picture, would appear very similar to what we see in pictures of 
Large Scale Surveys, such as the Lick or APM surveys ®. As he reduced the size of his 
gauge to smaller and smaller distances he would come to the conclusion that the Universe 
today is dominated by matter which, at large distances is in gravitational interaction. And 
that this matter seems to organize itself into bodies which roughly group themselves into 
solar systems, galaxies, groups of galaxies, and even larger structures. We believe we see 
the same as our tourist, and that we call the Cosmological Principle. 

At the larger scales he would probably recognize the Universe as a homogeneous object 
and therefore describe it with a Friedmann-Robert son- Walker metric. He would reproduce 
the many successes of cosmology based on this metric. As he went on to smaller scales 
he would find that (i) 'on all observable scales there are structures seen and significant 
anisotropics are detected'®. In fact, (ii) if he inferred the galaxy-to-galaxy correlation 
function (£,Gai( r )) from these surveys, he would discover^] that £(r) oc r~ 7 where § 7 ~ 
0(1.6 — 1.8), instead of 7 = 1 which is what one would naively expect^] for a homogeneous 
distribution of matter in a three dimensional space. 

Given these phenomenological facts, one may ask: (A) Is there a fundamental explana- 
tion for this power law behavior? (B) Can it be understood by using some basic scheme! 



The value of 7 seems to vary somewhat from survey to survey; for example, it is 1.8 for the Lick 
survey and about 1.6 for the APM. 

2 This can be understood as follows: as we will see below, the gas of galaxies can be put in a 
one-to-one correspondence with a 3-dimensional scalar field theory, and the galaxy-to-galaxy 
correlation function corresponds to the scalar field propagator. In units of length, a scalar field in 
a d-dimensional space-time has a canonical dimension of —(d/2 — 1). The above statement about 
7 follows at once. 
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This note aims at providing answers to the above questions within the framework of 
known physics. It rests upon the well-known observation H that the existence of a non- 
integer (anomalous) dimension (signalled by 7 7^ 1) is a tell-tale sign betraying the existence 
of both smaller length scales and fluctuations. 

We will apply the techniques of Statistical Mechanics to a system made up by many- 
galaxies (accounting for the 'smaller length scales' ) in gravitational interaction, and which 
are subject to fluctuations emerging from the intrinsic properties of the gravitational in- 
teraction of this many body system. The actual nature of the fluctuations needs not be 
specified here, but they could, for example, be related to the 'frictional' processes in grav- 
itational systems long ago considered by Chandrasekhar 11 or, as predicted by well known 
classical theorems of Poincare, to chaotic processes @ related to the many body nature of 
the gravitational system. 

We will obtain the partition function for this system and compute the two-point corre- 
lation function and its corrections due to the existence of fluctuations. 

Let us consider the continuous mass density, p(r), describing the spatial-distribution of 
galaxies. The deviation from an average density p is Sp(r) = p(r) — p. The galaxy-to-galaxy 
correlation function is defined @ as ^Ga/( r i — r j) = (^p( r i)^P( r j)) /p 2 > where the angular 
brackets mean that a suitable average has been taken. In the (justifiable) non-relativistic 
limit, the gravitational interaction energy for this system is then given by (assuming no 
expansion) 

H int = ^-G ( ( dhidhjp^i)-—^ — rpfo) . 

I * i 3 I 

As a first approximation, it is reasonable to consider the gas of galaxies as made up 
of discrete, spatially localized 'points' of (equal) mass mo, and we can set o the 'contrast' 
5p{Ti)/p equal to 1 if there is a galaxy at position and equal to -1 if there is a void. The 
interaction energy becomes 

Hi nt = —\Y1 m i \ r . _ r .| m j = ~\ 5 m i L iJ m j W 

ij I * j\ ij 



with the following natural correspondence (since we have assumed that all galaxies have a 
similar mass) between m i; a two- valued (±1) 'spin 1 variable s i: the density and the contrast: 

2 rrii/mo <-> p{vi)/p (2) 

Si^5p{Yi)/p. (3) 

Furthermore, mi and Sj are related by 

m i = m ^( y s i + l) (4) 

which happens to be analogous to the relationship between a lattice gas and an Ising magnet. 

This relationship, via a Hubbard-Stratonovich transformation 11 , allows one to map the 
'gas' of galaxies into a system described by a stochastic 1 -component (scalar) classical field 
0(r) in 3 dimensions, whose partition function Z[f3] can be readily calculated, and gives 

{m} 

= C J [d(j)) exp j-^ J (0(r) - H{r)) L~ x (r, r') (0(r) - H{r')) dvdr' + J dr log cosh [/?0(r)]| 

= J [dfle-M+M . (5) 

Here the function H(r) is given by H(r) = — 1/2/ dr'rriQL(r,r') and C is an inessential 
factor. The field <ft is the order parameter for this system. Although all direct reference to 
the original masses has disappeared, the physics described by the hamiltonian of Eq. (H) is 
completely equivalent to the original description. 

Because of Eq. ([!]) and (g), 



L-\r, r') = -{5{r - r')/{27rGm 2 )}V 

and 



2 

r ' 
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H[t/>, H] = ~\J (0(0 - H(r)) ^^V 2 (0(r) - (r)) dr 



dr log cosh [/30(r)] . 



27rGmg 

1 

As is well known @, the connected, two-point correlation function for the spin system 
and the field theory are the same. Furthermore, because of fluctuations in the field 0, its 
canonical dimension acquires an anomalous dimension and shifts away from its canonical 
value (cf. Footnote 2), such that when jr, — rj\ — > oo, the connected, 2-point correlation 
function for this hamiltonian, £(|r^ — r 3 -|), scales as @ o 

lim (s iSj )= lim ^(1^ - r^-l) ~ — 1 (6) 

|ri-rj|— >oo |r»— rj|-«>o \ r i ~ r j\ 

where the first equality follows from the equivalence between the 'spin' and ' field! descrip- 
tions of the system, d is the dimensionality of space (=3) and rj is the critical exponent for 
the pair correlation function, whose value (0.0198 - 0.064) (cf. the Table below) differs from 
zero due to the fluctuations in 4>(r). 
Because of Eq. (H) and @ 

tGal(\ri - Tj\) = (SiSj) (7) 

with the average computed using Eq. (|). 

Putting together Eqns. and (0), our calculation shows that for large separations, the 
galaxy-to-galaxy correlation function must scale as 



with @7 = d- 2 + r? between 1.0198 and 1.064. 

Thus far our calculation has been static, but the Universe is expanding, and the effects 
of expansion can (and do @) modify the values of critical exponents. For the correlation 
function, it is known from computer simulations in condensed matter physics § combined 
with dynamical scaling considerations, that time enters in the correlation function by altering 
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the argument of the correlation function from |r— to |rj— Tj\/L(t) where L(t) oc and ( is 
determined in computer simulations to be 1/3 for systems with a conserved order parameter, 
and 1/2 for systems with a non-conserved order parameter. Separation and time are related 
in an expanding Universe where, to a first approximation, in a matter dominated Universe 
the scale factor is proportional to t 2 ^ 3 . Putting together the expansion of the Universe and 
the dynamical critical phenomena effects as contained in (, the exponent in the galaxy-to- 
galaxy correlation function is modified from 7 = cf — 2 + 77 to (<i — 2 + 77) x (1 + 3(/2). That 
is, we finally get that the predicted (calculated) value for 7 will be[] between 1.530 and 1.596 
( = 1 Expanding) ^ we assume that the order parameter is conserved, and between 1.785 and 
1.862 (= ^Expanding) ^ we assume that the order parameter is not conserved. 



Method of Calculation 


IStatic 


"^Expanding 


1 Expanding 


Series Estimates 


1.056 ± 0.008 


1.584 ± 0.012 


1.848 ± 0.014 


0(e) 





1.5 


1.75 


0(6 2 ) 


1.0198 


1.530 


1.785 


0(e 3 ) 


1.037 


1.555 


1.815 


0(6 4 ) 


1.029 


1.543 


1.801 



These values are to be compared with the values inferred from the existing galaxy cata- 
logs, which range between 1.5 for the APM survey to 1.8 for the Lick survey B. 

Therefore, we see that the questions enumerated at the beginning of this note can find 
an answer within the framework outlined here. In addition, there is a clear and unam- 
biguous prediction: due to the Universal nature of the gravitational force, reflected in the 
interaction hamiltonian of Eq. (|]), the result we have obtained for the galaxy-to-galaxy 



5 Cf. Table. The static values for n are taken from Reference [10]. 



6 



correlation function must apply also to any other many-body-gravitational system, includ- 
ing the interstellar medium in our galaxy. This means that observations must confirm that 
the interstellar medium has a distribution whose correlation function scales with the same 
generic power law as galaxies. The only possible difference would be in the numerical value 
of the anomalous dimension, since for intergalactic gas clouds the size of the system is 
smaller, and therefore, renormalization group arguments tell us that the value of 7 interstellar 
is smaller than for systems of galaxies, where the coupling constant has had more distance 
to grow on its way into the IR fixed pointQ. 

Many questions remain. For example, is the order parameter conserved, or not? Can one 
use renormalization group techniques to also calculate the dynamical effect of the expansion 
of the Universe on 77, instead of appealing to phenomenological (albeit well substantiated) 
computer estimates to generalize the static values of 7 to dynamical values? Does the implied 
r-dependence of G(r) play a role in the physics of large scales? How do initial conditions 
impact on the correlation functions? How does one actually approach the disordered phase? 
These questions will be considered in a future paper. 



4 This follows by noticing that, because of conservation of probability (unitarity in field theory) 
n is positive; perturbation theory tells us that it is proportional to a power of Gm^. Also, in less 
than four dimensions, the latter coupling tends in the IR to the equivalent of the Wilson-Fisher 
fixed point, and away from the Gaussian UV-fixed point. 
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